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A recent description of an exact map for the equilibrium structure and thermodynamics 
of a quantum system onto a corresponding classical system is summarized. Approximate 
implementations are constructed by pinning exact limits (ideal gas, weak coupling), and 
illustrated by calculation of pair correlations for the uniform electron gas and shell struc¬ 
ture for harmonically confined charges. A wide range of temperatures and densities are 
addressed in each case. For the electron gas, comparisons are made to recent path inte¬ 
gral Monte Carlo simulations (PIMC) showing good agreement. Finally, the relevance for 
orbital free density functional theory for conditions of warm, dense matter is discussed 
briefly. 
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1. Introduction and Motivation 

A fundamental description of the thermodynamics (e.g., free energy) and structure (e.g., 
pair correlation function) for materials comprised of electrons and ions remains a challenge 
for many state conditions of current interest [lj. Typically the ions can be described by 
semi-classical methods due to their relatively large masses. In contrast, the electrons may 
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require an accurate description of strong quantum effects. Typical solid state conditions 
occur at temperatures well below the electron Fermi temperature for which the multi¬ 
plicity of zero temperature many-body theories and simulations are available. However, 
as the temperature is increased to several times the Fermi temperature such methods 
fail or become increasingly difficult to implement. At still higher temperatures, effective 
classical methods can be applied (e.g., molecular dynamics simulation [2], Monte Carlo 
integration, classical density functional theory [3], liquid state theory 0)- There is a long 
history of phenomenological attempts to extend these classical methods to lower temper¬ 
ature by including quantum effects in modified pair potentials j5j. More recently, these 
effective classical systems have been improved significantly by the inclusion of quantum 
effects in a modified classical temperature mm- The entire approach of constructing a 
classical system to replicate the thermodynamics and structure has been given a formally 
exact context from which more controlled approximations can be constructed PHO] . The 
objective here is to summarize briefly this latter work and to illustrate its utility by two 
applications: 1) the calculation of pair correlations for the uniform electron gas, and 2) 
the description of shell structure for charges in a harmonic trap. In both cases the em¬ 
phasis is on conditions ranging from classical to strongly quantum mechanical. The last 
section describes how this effective classical approach can be exploited to address current 
problems of ’’warm, dense matter” via orbital free density functional theory pQ. 

2. Definition of the Effective Classical System 

Consider a system of N particles in a volume V with pairwise interactions and an 
external single particle potential. The Hamiltonian is 

N 

H n = K + $ + J2v( q*), (1) 

i 

where K and $ are the total kinetic and potential energies, respectively. The form of 
the pair potential (f> (q,;, q^) and external potential v (q.;) is left general at this point. 
The equilibrium thermodynamics for this system in the Grand Canonical ensemble is 
determined from the grand potential S~2(/3 | /x, 0) 

I n, 0) = - In ^ Tr N e-^ K+ *~ S _ ( 2 ) 

N 

Here the local chemical potential /i(r) is defined by 

K t )=V-v(t), ( 3 ) 

and the operator n( r) representing the microscopic density is defined by 

N 

w(r) = ^5(r- qi ). (4) 

i =1 

The notation Q(P j /i, </>) indicates that it is a function of the inverse temperature /3 -1 = 
ksT and a functional of //(r) and the pair potential </> ( r,r'). 
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A corresponding classical system is defined with a classical grand potential 

m c (& I H', &) = - In X (5) 

N A ° J 

1 /2 

where A c = (2irf3 c h 2 /m) ~ is the thermal de Broglie wavelength and the inverse tem¬ 
perature of the quantum system 0 multiplies the classical grand potential. The classical 
grand potential is defined in terms of an effective inverse classical temperature 0 C , effective 
classical local chemical potential p c (r), and effective classical pair potential 0 C (r, r') . The 
classical system therefore has one undetermined scalar and two undetermined functions. 
These are defined by the following three conditions 


o /-o i .. j, \_n/a .. j.\ SQ C (P C \ p c ,0 c ) , _ 5Q(P | p,0) 

Qc\Pc | Pofyc) — Q\P I s I /3 C ,0c — 


H f c (r) 


<*M r ) 


( 6 ) 


sn c (f3 c I Pc, 0 C ) _ 5Q(P I p» . 

r 1 /_ A 1/3,M • 


<50 c (r,r') 


5<J)( r, r') 


( 7 ) 


An equivalent form for these conditions can be expressed in terms of the pressure, the 
local average density, and the pair correlation function 


Pc{Pc | Pc, </>c) = p{P \pA), n c (r;/3 c \ p c , <f> c ) = n( r;/3 | p, 0), 


( 8 ) 


p c (r, r'; /3 C \ p c , (j) c ) = p(r, r'; | p, 0). (9) 

In this way the classical system has the same thermodynamics and structure as that of 
the underlying quantum system. 

These definitions for /3 C , p c , and (f) c are only implicit and require inversion of the classi¬ 
cal expressions on the left sides of these equations to express them in terms of the given 
quantum variables /?, p, and <fi. Generally this is a difficult classical many-body problem. 
In addition, the inversion is expressed in terms of the corresponding quantum functions 
p(/3 | p, 0), n(r;/3 | p, 0), and g( r, r'; 0 j p, 0) which require solution to the original difficult 
quantum many-body problem. Hence it would appear that the introduction of a represen¬ 
tative classical system to calculate the thermodynamics of the quantum system is circular. 
However, it is expected that the inversion can be accomplished in some simple approxi¬ 
mation that incorporates relevant quantum effects and the resulting approximate classical 
parameters 0 c ,Pc and 0 C used in a more accurate theory or simulation to ”bootstrap” a 
better thermodynamics and structure. This is illustrated in the next two sections. 

3. Pair Correlations in the Uniform Electron Gas 

To illustrate the utility and effectiveness the effective classical system approach defined 
above, the calculation of pair correlations in the uniform electron gas is described in 
this section. The objective is to describe these correlations over the entire density and 
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temperature plane. In the classical domain this system is typically known as the one 
component plasma. The classical system is well-described by classical methods such as 
liquid state theory, classical Monte Carlo, and molecular dynamics simulation. Both fluid 
and solid equilibrium phases are now well characterized, including very strong coupling 
conditions. Consequently, there is a great potential to apply these approaches as well to 
quantum systems using the classical map. 

In this section the system of interest is the uniform electron gas at equilibrium. It is 
comprised of electrons in a uniform neutralizing background. The dimensionless temper¬ 
ature used here is the temperature relative to the Fermi temperature, t = 1/ [3eF, where 
the Fermi energy is defined by 6f = h 2 (37r 2 n) 2//3 /2m . Also the density dependence is 
characterized by the ratio of the average distance between particles relative to the Bohr 
radius, r s = r 0 /a,B, where Anrir q/3 = 1 and as = ti 2 /me 2 . The dimensionless space scale 
is r*= r/r 0 . The classical pair correlation function at uniform equilibrium depends only 
on the relative coordinate so 


g c (r,r']P c \ // C ,0 C ) = g c {\r* - r*'| ,r s | 0*), (10) 

where 0* = /3(f)* = r/|r* — r*'|. Here T = f3e 2 /rQ is the classical Coulomb coupling 
constant. In terms of t,r s it is T = 2 (|-7r) 2 ^ 3 r s t _1 . Note that the functional g c (r,r s j •) 
is independent of t. In contrast, the quantum pair correlation functional depends on both 
r s and t 

g{ r, r '; P I he, 0 C ) = g{\r* - r*'| , r s , t \ <j>*). (11) 

Now, using the equivalence of the classical and quantum pair correlation functions (|9]) 
the classical functional can be inverted to give 0* 

0c (r*,r 8 ,t) = g~\r*,r s \ g). (12) 

This is the formally exact definition of the effective classical pair potential. The prac¬ 
tical procedure is to evaluate this in some reasonable, simple approximation and then 
’’bootstrap” the result in a more sophisticated approximation to g c (r*,r s \ <fr*). First, it is 
required that the limit of non-interacting particles be given correctly, so the potential is 
written as 


0c ( r *,r s ,t) = 0* (o) (r*,r s ,t) + A(r*,r s ,t). (13) 

Here 0^°'* (r*,r s ,t) is an effective pair interaction chosen such that its classical pair cor¬ 
relation function is the same as that for the quantum system with no Coulomb interac¬ 
tions, g^ (r*,r s ,t) . The second term, A (r*,r s ,t), replaces the Coulomb interaction by 
a corresponding classical pair interaction incorporating the quantum effects. Here, it is 
constrained to be exact in the weak coupling limit. Classically, the latter corresponds to 
the potential becoming the same as the direct correlation function 


0* (r*,r s ,t) ->• —c (r*, r s , t ), 0* (o) (r*,r s ,t) ->• -c (0) (r*,r s ,t) 


(14) 


where the direct correlation function is defined by the Ornstein-Zernicke equation 


c (r*,r s , t) — g (r*, r s ,t) — 1 — — / dr'c (r*', r 3 , t) (g (|r* « r*'| , r 8 , t) - 1). (15) 
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The quantum weak coupling limit is the random phase approximation. Therefore © is 
calculated by inserting the finite temperature random phase approximation, g RPA ( r*,r s , t ), 
on the right side. The resulting approximate effective pair potential is now (13) with 


(f)* c (r*,r s ,t) ~ (r*,r a ,t) - c HFA (r*,r a ,t) + c w (r*,r a ,t) . 


RPA 


.(0) 


(16) 


Clearly this incorporates the ideal gas and weak coupling limits, without being restricted 
to either. 

The qualitative differences of this effective classical potential from the underlying Coulomb 
potential of the quantum system are two fold. First, the divergence at r* = 0 is removed, 
i.e. (j )* (r* = 0,r s ,t) is finite. Second, for large r* the potential is also of the Coulomb 
form, but with a different amplitude 


(p* c ( r *,r a ,t) —>r e (t,r s ) r* \ (17) 

The classical Coulomb coupling constant F(t,r s ) has been replaced by the effective 
quantum coupling constant T c (t, r s ) 


T e (t,r a ) 


2 

PUujp cotli (/3hu p / 2) 


r, 


T (t , r a ) = /3q 2 /r 0 



(18) 


Here /3hu p = (3h (Airnq 2 /m) 1 ^ =4 (2v / 3vr _2 ) 1 ^ 3 ry 2 /3f is the dimensionless plasma fre¬ 
quency. 

With the pair potential determined in this way the pair correlation function g (r*,r s ,t) 
can be calculated beyond the ideal gas and weak coupling conditions using, for example, 
molecular dynamics or classical Monte Carlo simulation. Here the results are illustrated 
using an integral equation from liquid state theory or classical density functional theory. 
It is the hypernetted chain approximation (HNC) 


In g{r*,r a ,t) = -</>* (r*,r a ,t) - c(r*,r a ,t ) + (g(r*,r a ,t) - 1) (19) 

This equation together with the Ornstein-Zernicke equation ( ]l5| ) provides a coupled 
set of equations for both c(r*,r s ,t ) and g(r*,r s ,t ). Figures 1 and 2 show the results in 
comparison with recent path integral Monte Carlo (PIMC) simulations [TT] at r s — 1 and 
6 for a wide range of t. Clearly there is quite good agreement with this benchmark data 
using this standard liquid state classical theory modified only by the quantum effects in 
the modified pair potential. For additional details and other values for r s , t see reference 

ra- 


4. Charges in a Harmonic Trap 

For a second application of the classical map consider N charges in a harmonic trap. 
The classical HNC of the last section, extended to this inhomogeneous system [13], has 
been shown to give an accurate description of the radial density profile for this system 
ni- Of particular interest is the quantitative description of shell structure that occurs for 
classical strong coupling conditions. In this section, that approach is extended to include 
quantum effects. 
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Figure 1. Pair correlation function g(r*) for 
r s — 1 at t — 0.5, 1, and 8. 



Figure 2. Pair correlation function g(r*) 
for r s = 6 at t = 0.5, 1, and 8. 


The HNC theory for the density profile together with Q and (J9]) leads to 


111 (n (r) Aj) = - ffv c (r) + / dr'c(|r —:r'| ,/i c ,ft)n(r'). 


( 20 ) 


Here c(r, /x c , f3 c ) is the same direct correlation function as described in the previous section. 
Also, v c (r) is the effective classical trap which includes quantum behavior of the given 


harmonic trap. The approach is to determine u c (r) approximately by inverting (20) for a 
given approximate quantum density n (r). Two approximations are compared here. 

The first approximation is to invert (20) for non-interacting particles in a trap. The 


quantum effects in this case are entirely due to exchange symmetry 


do) 


(r*) = c — In (n**- 0 ) (r*)) + / dr* 7 c^(|r* — r* 7 | ,r s ,t)n*^ (r* 7 ) 


( 21 ) 


Here u*(r*) = Ar c ( r ) and the superscript 0 on a property denotes its ideal gas value. Also, 
c is a constant that only sets the normalization of the density profile. The calculation 
of n*^ (r*) is straightforward in terms of the harmonic oscillator eigenfunctions, but 
for the case considered here (N = 100) it is found that the local density approximation 
(finite temperature Thomas-Fermi) is quite accurate. An important qualitative feature 
of n (r*) is its vanishing at a finite r* as t —» 0. This leads to the formation of a 
hard wall in the effective classical trap. It is well known that such hard walls produce 
shell structure in classical mechanics, so this represents a quantum origin for new shell 
structure independent of Coulomb correlations. 


The second approximation is to invert (20) with mean field quantum Coulomb correla¬ 
tions 


C {H) ( r*) = c-ln 


n 


.(H) 


(r*)) + J dr*' (V°)(|r*—r* 7 | ,r s ,t) - n*™ (r *'). (22) 
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The density n*^ (r*) is calculated from quantum density functional theory without ex¬ 
change or correlation (Hartree approximation), and again using the local density approx¬ 
imation. It gives a qualitative change from the ideal gas form (21) since the system is 
considerably expanded by the Coulomb repulsion. The hard wall is mitigated and result¬ 
ing effective classical trap potential has a more harmonic form. 

Figure 3 shows the density profiles obtained from (20) using (21) or (22). It is seen 
that the strong shell structure from the ideal gas hard wall is removed when Coulomb 
interactions are included. Also shown on this figure are the effective classical trap poten¬ 
tials for the two cases. The result from (21) shows a kink which is a precursor of the hard 
wall at t — 0. Although this deviation from harmonic is small, it is sufficient to generate 
a large shell. In contrast, the result from (22) is more nearly harmonic and has only 
the shell due to Coulomb correlations already present in the purely classical calculation 
(no quantum effects) [13]. The lesson from this comparison is that quantum effects on 
the effective classical trap potential determined without Coulomb interactions lead to a 
false mechanism for shell structure. The more realistic mean field quantum determination 
does not have this shell structure and provides a quite different density profile. A more 
complete discussion of this comparison and results for a wide range of r s ,t will be given 
elsewhere 



Figure 3. Density profiles (thick lines) and 
effective classical trap potentials (thin lines) 
using the Hartree and ideal gas approxima¬ 
tions. 


5. Discussion 

The use of an effective classical system to describe quantum effects has been shown 
to provide a practical tool by the two examples of the previous sections. Although the 
results are quite good for the methods used here to determine the effective pair potential 
and trap potential, and the HNC implementation of the classical statistical mechanics, 
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improvements in both remain to be explored. For example, the limitations of the HNC 
theory can be eliminated using these same potentials in molecular dynamics or Monte 
Carlo simulations. 

The comparisons of the last section show that incorrect results can be obtained if the 
input quantum mechanics for the effective potentials is not sufficiently representative of 
the real system. For the electron gas, the effective potential has a Coulomb tail whose 
amplitude is constrained to satisfy an exact sum rule. It would be useful to have exact 
constraints for other systems as well, to assure applicability over a wide range of the 
parameter space. 

As noted in the introduction, there is strong current interest in systems of electrons 
and ions to describe conditions of warm, dense matter [1|. Such systems are described 
by molecular dynamics simulation of the ions whose forces are calculated from a den¬ 
sity functional theory (DFT) for the electrons at each time step. Traditionally, the DFT 
calculation is performed within the Kohn-Sham approach requiring a self-consistent di- 
agonalization of an effective single electron Hamiltonian to construct the local density. 
At temperatures approaching the Fermi temperature, the number of relevant states (or¬ 
bitals) becomes large and the calculations are no longer practical. A resolution of this 
problem is to forgo the Kohn-Sham method and return to the original form of DFT with 
a single Euler equation for the local density determined from the free energy as a known 
functional of the density. The primary difficulty is hireling the non-interacting free energy 
as a functional of the density, which remains an unsolved problem in the quantum theory. 
However, its classical counterpart does not have this difficulty — the non-interacting free 
energy is known as an explicit functional of the density. Hence, an implementation of 
the effective classical system as described here, together with classical DFT, provides the 
desired orbital free DFT. 

To see how this might be implemented consider a system of N e electrons and N t positive 
ions with charges Z and positions {R a }. For charge neutrality N t Z = N e . This can be 
viewed as an electron system in the external potential of the ions 


A -Ze 2 

u ( r ) = Z. EFT 

a=l ' 


R„ 


(23) 


Return to 


for the corresponding local electron density, where now u c (r) is the effective 


classical potential corresponding to (23). A reasonable, realistic determination of u c (r) 
might be given by (22) in the form 


v< HF \r* 


= c — In (n« HF > (r*)) 


dr 


*' l c (o) 


(l r * 


,r s ,t) 


r (t,r a ) 


n 


*{HF ) / */ 


(O. 


(24) 


where now n*^ HF > (r*) is the Hartree-Fock electron density for the given array of electrons. 
While determination of n*^ HF ^ (r*) is still non-trivial it is a practical problem, and then 
use of v*^ HF \r*) in (2o| gives the desired orbital free DFT for the electrons. 

Dharma-wardana has proposed a more complete application of the classical DFT for 
both the electrons and ions [7j, eliminating the molecular dynamics simulation for the 
ions. An additional effective classical electron - ion potential must be determined in this 


case. 
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